Quantization of Superflow Circulation and Magnetic Flux with a Tunable Offset 
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Introduction: In superfluid '*He, as well as in conventional 
superconductors, for any configuration of the appropriate 
order parameter, and for any closed path through nonsin- 
gular regions of it, the superflow circulation k adopts one 
of only a discrete set of values. Quantum-mechanical in 
origin , this remarkable phenomenon is referred to as 
the quantization of circulation. In "^He, circulation occurs 
at integral multiples of h/M, M being the mass of the 
''He atom. (In conventional superconductors, M is the 
Cooper-pair mass). The magnetic flux enclosed in a ring 
of conventional superconductor is also quantized, in units 
of $0 (= 2.07 X 10^^^ We). In superfluid ^He, however, 
K is not necessarily quantized, owing to the anisotropy of 
certain of its superfluid phases ||^,02p|. 

For •^He-A, the standard view ^[ is that despite the 
absence of circulation quantization in the hulk, integral 
quantization, precisely as in ^He, occurs at surfaces, ow- 
ing to the anchoring of the relative orbital angular mo- 
mentum vector 1. In fact, we shall see that quantization 
at surfaces is by no means the rule, instead being the 
exception, with integral quantization being yet more ex- 
ceptional. Indeed, the aim of this Letter is to demon- 
strate that, although quantization is not generic, there 
do exist experimentally achievable order parameter con- 
figurations for which circulation is quantized, but not to 
conventional (i.e., integral) values. Instead, what occurs 
is offset circulation quantization, i.e., integrally- spaced 
hut nonintegral circulations: ■ ■ —1+10,0 + uj,l +10,- ■ •, 
for all paths. The offset lo can be continuously tuned 
by varying the sample shape or certain external fields, 
both of which can determine the texture of the order 
parameter. Moreover, as we shall see, the Aharonov- 
Bohm-like physics responsible for offset circulation- 
quantization in ■^He-A also leads to the possibility that 
a nonintegral numher of magnetic flux quanta could 
be enclosed in a ring of unconventional superconduc- 
tor. Furthermore, if the order parameter is in equilib- 
rium then offset quantization has the striking corollary 



that arbitrary nonintegral equilibrium circulations can 
be obtained, in contrast with the standard possibilities: 
0,1/2,1 (bulk); 0,1/2,1,3/2,... (surface); see, e.g., ||. 

We proceed as follows. First, we establish a gen- 
eral relationship between circulation, parallel transport 
and order parameter anholonomy. Next, we introduce 
Aharonov-Bohm-type concepts to this area of physics. 
Finally, we discuss experimental configurations for su- 
perfluid ■^He-A, and propose tests involving the offset of 
magnetic-flux-quantization in unconventional supercon- 
ductors and the phase-shift of Little-Parks oscillations 

i- 

Geometric phase: The origin of the phenomena consid- 
ered here is the geometric phase [^,0, i.e., the phase ac- 
quired by Cooper pairs as they propagate through the 
condensate that they self-consistently form, their rela- 
tive (orbital or spin) angular momentum following the 
local orientation determined by the condensate. This ge- 
ometric phase is an "anholonomy: the geometrical phe- 
nomenon in which nonintegrability causes variables to 
fail to return to their original values, when others, which 
drive them, are altered round a cycle" ||l^. Indeed, as 
we shall show, the circulation computed around any path 
in an anisotropic superfluid is a natural anholonomy in 
the parallel transport of the triads that characterize its 
order parameter. 

Consider superfluid ^He-A, which is characterized by 
an order parameter matrix d^j of the form: 



Ao d^[mj 



exp IX. 



(1) 



Here, Aq is the (real) gap parameter for the conden- 
sate, dfj, is a (real) unit-vector associated with the spin 
angular momentum sector, {1, m, n} is an orthonormal, 
right-handed, triad of (real) vectors associated with the 
relative orbital angular momentum sector, and x is a 
(real) phase. This parametrization of d^j in terms of 
X and {1, m, n} is not unique: d^j is invariant under 
the simultaneous operations of rephasing x (by 5x) and 
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rotating {m, n} about 1 (by —Sx)- Due to the bro- 
ken relative gauge-orbital symmetry (for reviews, see 
Refs. the superfluidity is manifested not by 

a complex scalar, as in "^He, but instead by the complex 
vector t = (m -I- m) /V2. 

This recognition, that the parametrization oidnj in 
terms of x and {1, m, n} is not unique, is crucial At 
any point in (real) space one may choose, e.g., x = 0, and 
thus select a particular {1, m, n}. A choice for x and the 
pair {m, n} at some other point may, if one wishes, then 
be determined as follows: given a path between the two 
points, the orientation of {m, n} at the second point is 
obtained by the parallel transport (defined below) of the 
pair over the (curved) surface of orientations of 1 that are 
encountered along the real-space path. However, owing 
to the nonintegrability of this parallel-transport law, the 
orientation of {m, n} cannot be expressed as a single- 
valued function of 1 over the entire 1-spherc. Thus, if the 
(real-space) path is closed then, although 1 returns to its 
original value, m and n and therefore x will, in general, 
not (i.e., there can be anholonomy). 

The parallel transport of triads over the 1-sphere is 
accomplished via the connection 

Imt*Vitj=0, (2) 

the line integral of which turns out to be the solid angle 
swept out by 1 as the path is traversed , and gives the 
anholonomy mentioned above. 

Circulation and geometry: Having discussed geometric 
issues associated with the superfluid ^He order pa- 
rameter matrix, we now turn to the relationship between 
these issues and the circulation k. The local superfluid 
velocity v^^' fll] can be expressed in terms of the orbital 
triad {1, m, n} (or, equivalently, the complex vector t) 
via 

2Mv(^' /h = rrijVr rij + Vr X = Im Vr tj + Vr X- (3) 

That v^''-' can be interpreted as the superfluid velocity 
follows from its behavior under galilean transformations: 
if t ^ t' = texp {2iMSv ■ r/h) then v(") ^ + Sv. 
A local rephasing of x and rotating of {m, n} (such that 
dfij is unchanged) amounts to a gauge transformation, 
and thus leaves v'^*') unchanged. The velocity v^**) has 
the form of a gauge potential, like the London current in 
a superconductor [1^]. In the present context, this gauge 
potential is referred to as the Berry connection. As d 
is a real vector in the A-phase of '^He, it produces no 
contribution to v^*^^ , even if d is inhomogeneous. 

Now, K is the line integral around a closed path Cr in 
real space: k = <f^ dr ■ v^**) . By using Eq. (^) k can be 
expressed as 

K = {n/2M) i dr ■ {lm<* Vr tj + Vr x} (4a) 
= {n/2M) i d\ ■ {Im t* Vi tj+Vix), (4b) 



where C/ is the contour on the 1-sphere traced out by 
the orbital vector 1 as the path Cr is circumnavigated. 
By adopting the parallel-transport connection, Eq. (||), 
we see that the penultimate term in Eq. (^b|) vanishes, 
so that K is determined by the anholonomy x(C/) in Xj 
which is necessary to compensate for the anholonomy in 
{m, n} induced by the parallel-transport connection: 

ti={n/2M)x{Ci). (5) 

Hence, we see that the superfluid circulation is given by 
the anholonomy in the parallel transport of the order pa- 
rameter, and that k is the flux associated with the Berry 
connection. The computation of x(C/) fl^ gives 

x{Ci) ^ niCi) +2np, (6) 

where ^l{Ci) is the solid angle subtended by Ci at the 
center of the 1-sphere and p is an integer. Although ^l{Ci ) 
is deflned only modulo 47r, the term 27rp (including odd 
p) results in the standard circulation quantum. 

For a nonsingular order parameter matrix in a simply- 
connected sample, Eq. (^) is of course the well-known Ho 
circulation theorem, obtained by Ho in a different manner 
p6| . The parallel-transport approach adopted here yields 
an extension to multiply-connected samples, flrst conjec- 
tured in Ref. |0. The extra term, 2p7r, allows for sin- 
gularities, which render the sample multiply-connected. 
We emphasize the anholonomic essence of the circulation, 
even in simply-connected samples. We note that Eq. m) 
holds whether or not the distributed vorticity V x v*^' 
of the superflow, which is given by the Mermin-Ho p^ ] 
equation Vr xv'*'' = {h/AM)eijkh (Vr x Vr Zfe) , is zero. 

We now come to our primary observation. In a 
multiply-connected sample from which distributed vor- 
ticity is absent the circulation k is quantized with an off- 
set, in striking contrast with superfluid *He. As we elu- 
cidate below, this amounts to Aharonov-Bohm physics 
stemming from geometry. The offset is determined by 
the solid angle fl{Ci) via Eqs. (||) and (^. Further- 
more, the offset of k results in truly stable (i.e., equi- 
librium) persistent flow. For a system with radius of 
order 1 mm the speed of this equilibrium flow will be of 
order 3 x 10~'^mms~^. It is desirable to exclude dis- 
tributed vorticity because whenever the area between 
nearby (real-space) paths is penetrated by distributed 
vorticity, a continuous variation of a path will lead to 
a continuous variation of k. (The presence of a small 
amount of distributed vorticity will smear the sharpness 
of the quantization, but not destroy it altogether.) By 
contrast, if the sample is simply-connected and there is 
distributed vorticity then k is not quantized (i.e., a con- 
tinuous variation in the path leads to a continuous varia- 
tion in k). However, contrary to the standard view, quan- 
tization is absent not only in the bulk but also at surfaces, 
owing to their geometry. For example, for surfaces with 
nonzero gaussian curvature (such as hyperboloids) there 
is no quantization. (If the sample is simply-connected 
and there is no distributed vorticity then k is zero.) 
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What is the connection between Aharonov-Bohm 
physics and offset circulation-quantization? Due to the 
structure of the ^He-A order parameter, the distributed 
vorticity of the superflow velocity field v'^*') is, in gen- 
eral, nonzero. Hence, in a generic sample a contin- 
uum of values of the circulation can be found. If the 
distributed vorticity vanishes and the sample is simply- 
connected then for all circulation-paths the circulation 
is zero. However, if the sample is multiply-connected 
then, even if the distributed vorticity vanishes, the circu- 
lation need not be zero. This is Aharonov-Bohm physics 
pgf . A close analogy with Aharonov-Bohm-type phe- 
nomena from mesoscopic physics is evident. Consider a 
thin- walled normal metal cylinder in a homogeneous ax- 
ial magnetic field . If the wall is sufficiently thin that 
a negligible amount of magnetic flux penetrates the metal 
itself then the conductance of the cylinder oscillates with 
the magnetic flux threading the interior of the cylin- 
der because all topologically-equivalent Feynman (elec- 
tron) paths enclose identical magnetic flux. If the wall 
is thicker, so that a non- negligible amount of flux pene- 
trates the metal, then the amplitude of the oscillations 
is diminished, as topologically-equivalent paths now en- 
close differing fluxes. In '^He-A, the superfluid velocity 
plays the role of the magnetic vector potential, and the 
distributed vorticity of the superflow plays the role of the 
magnetic field. A quantization offset, the analogue of a 
threading flux, requires a nontrivial ^{Ci), which can be 
maintained, e.g., by a surface or an external field that 
induces an appropriate texture in 1. 
Illustrative settings: We now illustrate the general ideas 
presented above by discussing two novel and experimen- 
tally feasible settings in which offsets may be observed. 
We note that several authors have considered a variety of 
textures for '^He-A, with and without singularities; see, 
e.g., Refs. |,|l|,|l]j22|j2|. (Equilibrium imphcations of 



geometric phases in the mesoscopic context have been 
discussed in Ref. [^.) In the first setting, a sample of 
■^He-A is contained in the cavity between a pair of long, 
coaxial, truncated cones of vertical angle a. The texture 
is maintained by surface orientation effects: the vector 
1 tends to be oriented perpendicular to walls ||2^]. We 
suppose that the separation between the cones is suffi- 
ciently small that the texture adopts the form shown 
in Fig. |l| and defined in the caption. Then, from Eq. (||) 
we see that 



{fi/2M p)e^ sin a, 



(7) 



and V X v^''^ = 0. Moreover, for this conical texture 
V X 1 = (see fljl). The offset of k for any path sur- 
rounding the axis once is given by {h/2M)il. (C/), where 
(Ci) is the solid angle subtended by the normal to the 
conical surfaces. This offset induces a persistent equilib- 
rium current. 

A persistent ground-state current in a conical container 
can be understood in terms of angular momentum con- 
servation. To see this, consider a cylinder (a — 0). Then 



ri (C;) = 2tt, the axial component of net relative orbital 
angular momentum vanishes, and k = 0. Now imagine 
varying a so that the cylinder becomes a cone, all the 
while maintaining cylindrical symmetry. Then, no axial 
torque acts during this variation, and the increase in the 
axial component of the net relative orbital angular mo- 
mentum is compensated by the appearance of center-of- 
mass orbital angular momentum and, hence, circulation. 

We remark that although the truncated cone and the 
cylinder are topologically equivalent they are, of course, 
geometrically inequivalent, this inequivalence being char- 
acterized by the vertical angle a. We hope that the com- 
plication arising from the truncations of the cones will, 
at least in the case of long cones, not be too severe. 

The second setting concerns a magnetic-field-induced 
spin-texture in the Ai phase of superfiuid '^He. As this 
phase is spin-polarized, so that only spin-projection +1 
Cooper pairs form the condensate, mass supercurrents 
imply spin supercurrents. The order parameter matrix 
has the form df^j — ^A|| (d^ -I- ie^) (rrij + irij^e^^. Here, 
is the (real) gap parameter, {f , d, e} is an orthonor- 
mal triad of vectors associated with the spin angular mo- 
mentum sector, and x is a phase. In the Ai phase, the 
orbital and spin sectors are identical in structure, and 
the superfluid velocity is given by 

vT = (;i/2M){lm (t* Vr t, + T* Vr T,) + Vr x}, (8) 

where t = {d + ie) / \f2. The result for k then reads: 



(9) 



K = (n/2M){l7i(C;) +%(C/) +27rp}. 



Both 1- and f-textures may lead to the offset of k and 
to a persistent equilibrium current. Consider a sample 
contained in the cavity between a pair of coaxial cylin- 
ders. A static, uniform magnetic field Bm is applied to 
stabilize the >li-phase (see [^||J|]). In this phase, the 
spin-orbit interaction ||^ tends to align (or anti-align) f 
and 1. If 1 is strongly anchored by the surface, as it would 
be in a sufficiently narrow cavity, then for Bm = B^Bz we 
have f = Bp cos P — e^sin/3, where tan/? = B^/Bgo and 
i3so(~ 2.8 mT) is the magnitude of the effective magnetic 
field ^ due to the spin-orbit interaction The sohd 
angle subtended by f is then given hynf = 27r(l-sin/3). 
This solid angle and, correspondingly, the persistent cur- 
rent, can be tuned by varying of B^. By tilting Bm away 
from the cylinder axis, flf is reduced. In the present set- 
ting, the sense of the current is determined by the mag- 
netic field, which explicitly breaks time-reversal symme- 
try. In the first setting, by contrast, time-reversal sym- 
metry is spontaneously broken (i.e., 1 may point inward 
or outward), the sense of the current being determined 
by the exhibited choice. Spin textures may also be con- 
trolled via the interplay between Bm and a second field 
(produced, e.g., by a axial current-carrying wire) that is 
inhomogeneous on the length-scale of the sample. 
Flux- quantization: We now turn to implications of geo- 
metric phases in the context of flux-quantization in su- 
perconductors with nonscalar order parameters. In a 
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superconductor that exhibits a spin-triplet pairing state 
IPSj and a corresponding order parameter matrix A^j the 
current is given by 

j = {2eh/M) Im A;^- + 2eA/h) A^^- , (10) 

where e is the electronic charge. Following Ref. we 
consider a ring of this superconductor, integrate over a 
path deep inside the ring (where the current is zero), 
and obtain the result that the enclosed magnetic flux 
is quantized with an offset. Similarly, we anticipate a 
phase-shift in Little-Parks oscillations Moreover, it 
would be interesting to test these ideas in the context 
of high-temperature superconductivity. Related effects 
may occur in the context of neutron star physics (see, 
e.g., Ref. [§, sec. 6.2.5). 

Conclusions: We have considered implications of geo- 
metric phases in the context of nonscalar superfluid- 
ity and superconductivity. We have shown that the 
geometry of the order parameter can have fascinat- 
ing quantum-mechanical ramifications, especially in the 
realm of circulation- and flux-quantization. In particu- 
lar, we have shown that, in the absence of distributed 
vorticity, the possible values of the circulation are deter- 
mined solely by the anholonomy in the parallel transport 
of the order-parameter triad along the circulation-path. 
We have also discussed certain experimental settings in 
which the effects proposed here might be observed. Thus, 
we have demonstrated that, even when irrotational, su- 
perflow in '^He-A is strikingly different from that in "^He 
and "^He-B. The present work may have implications for 
the experiments reported in Ref. pg| ] 
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FIG. 1. a) Concentric truncated cones of vertical angle a 
with superfluid "^He-A in the cavity. The texture is given 
by X — and (for cyhndrical polar coordinates {p, ip, z}) 
{1, m, n} = {— Gp cos a + ez sin a, Bp sin a + Bz cos a, e^}. The 
texture has 1 pointing from outer to inner cone [M. C,. is a 
path along which k is computed, b) 1-sphere with tangent 
plane containing {m, n}, and Ci is the image of Cr- 
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